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ABSTRACT 

We numerically investigate the migration of dust particles with initial orbits close to those of 
the numbered asteroids, observed trans-Neptunian objects, and Comet Encke. The fraction 
of silicate asteroidal particles that collided with the Earth during their lifetime varied from 
1.1% for 100 micron particles to 0.008% for 1 micron particles. Almost all asteroidal particles 
with diameter (i>4 microns collided with the Sun. The peaks in the migrating asteroidal dust 
particles' semi-major axis distribution at the n:(n+l) resonances with Earth and Venus and 
the gaps associated with the 1:1 resonances with these planets are more pronounced for larger 
particles. The probability of collisions of cometary particles with the Earth is smaller than for 
asteroidal particles, and this difference is greater for larger particles. 



INTRODUCTION 

The orbital evolution of interplanetary dust particles originating from the trans-Neptunian 
belt [l]-[5], Halley- and Encke-type comets [6], [7], and asteroid families and short-period comets 
[8]-[10] has been numerically simulated by several authors in order to understand the distribution 
of dust in the solar system. Observationally, the infrared satellites IRAS and COBE detected 
thermal emission from dust outside the Earth's orbit [11]-[13] while Pioneers 10 and 11 measured 
the interplanetary dust population directly using impact detectors [14]. The main contribution to 
the zodiacal light is from particles with diameters of about 20 to 200 fim. It has been estimated 
that 20,000 to 40,000 tons of micrometer- to millimeter-sized dust particles fall to the Earth 
every year with the mass distribution of dust particles peaking at about 200 ^m in diameter 
[15]. Liou et al. [7] showed that the observed shape of the zodiacal cloud can be accounted 
for by a combination of about 1/4 to 1/3 asteroidal dust and about 3/4 to 2/3 cometary dust. 
The hypothesis that the 'kuiperoidal' dust (i.e., the dust formed in the Edgeworth-Kuiper belt) 
is dominant in the outer solar system explains the Pioneer and Voyager data fairly well, while 
fitting a sample of the COBE/DIRBE data indicates that the kuiperoidal dust contributes as 
much as 1/3 of the total number density near the Earth (for a review, see [4]). Ozernoy [4] 
considered 1 m and 5 m particles while constructing the brightness of a disk of asteroidal, 
cometary, and kuiperoidal dust particles which fit the COBE/DIRBE data. 
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Analysis of the Pioneer 10 and 11 mctcoroid detector data [16], [17] showed that a population 
of 10^^ and 10^^ g (~10 /im) particles has a constant spatial density between 3 and 18 AU. 
Dust grains released by Halley-type comets cannot account for this observed distribution [6] , but 
trans- Neptunian dust particles can [2]. Using trans- Neptunian dust particles, Moro-Martin and 
Malhotra [5, Fig. 13] found that the surface (not spatial) density was approximately constant 
in this region. 

Liou et al. [1] noted that interstellar dust particles with an average size of 1.2 //m can destroy 
dust particles formed in the solar system and that the collisional lifetimes for 1, 2, 4, 9, 23 /im 
particles are 104, 49, 19, 4.8, 0.86 Myr, respectively. In these size ranges mutual collisions are 
not as important as collisions with interstellar grains. Moro-Martin and Malhotra [5] concluded 
that collisional destruction is most important for kuiperoidal grains between 6 //m (9 //m in 
Liou et al.) and 50 jj.ni. Particles larger than 50 //m may survive because interstellar grains 
are too small to destroy them in a single impact. For silicate particles 1-40 /um in diameter, 
the sublimation temperature (^^1500 K) is reached at i?<0.5 AU, but for water ice particles the 
sublimation temperature (-^100 K) is reached at 27, 19, 14, 10, and 4.3 AU for the sizes of 3, 6, 
11, 23, and 120 //m, respectively [1],[5]. 

In the present paper we consider a different set of initial orbits for the dust particles than 
the aforementioned authors and, for the first time, investigate the collisional probabilities of 
migrating particles with the planets based on a set of orbital elements during their evolution. 
We also present plots of the orbital elements of the migrating particles. 

MODELS FOR ASTEROIDAL DUST PARTICLES 

Using the Bulirsh-Stoer method of integration, we investigated the migration of asteroidal 
dust particles under the influence of planetary gravity (excluding Pluto), radiation pressure, 
Poynting-Robertson drag, and solar wind drag for values of the ratio between the radiation 
pressure force and the gravitational force /3 equal to 0.004, 0.01, 0.05, 0.1, 0.25, and 0.4. Burns 
et al. [18] obtained P=0.573Qpr/{ps), where p is the particle's density in grams per cubic 
centimeter, s is its radius in micrometers, and Qpr is the radiation pressure coefficient {Qpr 
is close to unity for particles larger than 1 fim). For silicates, the f3 values 0.004, 0.01, 0.05, 
0.1, 0.25, and 0.4 correspond to particle diameters of about 100, 40, 9, 4, 1.6, and 1 microns, 
respectively. Silicate particles with /3 values of 0.01 and 0.05 have masses of 10~^ g and 10~^ 
g. For water ice, our /3 values correspond to particle diameters of 300, 120, 23, 11, 6, and 3 //m 
[18]. We assume the ratio of solar wind drag to Poynting-Robertson drag to be 0.35 [5], [6]. 
The relative error per integration step was taken to be less than 10~*. 

The initial positions and velocities of the asteroidal particles were the same as those of the 
first numbered main-belt asteroids (JDT 2452500.5), i.e., dust particles are assumed to leave the 
asteroids with zero relative velocity. We considered A''=500 particles for each /3>0.05, A^=250 
for (3=0.01, and iV=100 for /3=0.004. With /3>0.01 in each run we took iV=250, because for 
A''>500 the computer time per calculation for one particle was several times greater than for 
N=250. The simulations continued until all of the particles either collided with the Sun or 
reached 2000 AU from the Sun. The lifetimes of all considered asteroidal dust particles were 
less than 0.8 Myr, except for one particle with a lifetime of 19 Myr. We also made similar runs 
without planets to investigate the role of planets in interplanetary dust migration. 

COLLISIONS of ASTEROIDAL PARTICLES WITH PLANETS AND THE SUN 

In our runs, planets were considered as material points, but using orbital elements obtained 
with a step dt of <20 yr (dt=10 yr for /3 equal to 0.1 and 0.25, and dt=20 yr for other values of 
/3), similar to [19] we calculated the mean probability P=Pj:/N (Ps is the probability for all N 
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Table 1. Values of T, Tj, T^**^, r^«^, T^^^"^, T^^g (in Kyr), P^, and Psun obtained for asteroidal 
dust particles for several values of /? (Venus=V, Earth=E, Mars=M) 







V 


V 


E 


E 


M 


M 














Psun 


Pr 


T 


Pr 


T 


Pr 


T 


Tj 




r-pmax 


rpmin 
-'2000 


-^2000 


0.004 


1.000 


12783 


40.5 


11350 


90. 


1204 


220. 





348 


932 






0.01 


1.000 


1534 


19.2 


1746 


44.2 


127 


99.9 





142 


422 






0.05 


0.996 


195 


4.0 


190 


8.1 


36.7 


20.5 





30 


89 






0.1 


0.990 


141 


2.4 


132 


4.8 


16.4 


12.0 


2.21 


16 


44 


138 


793 


0.1* 


0.990 


366 


2.4 


279 


4.8 


20.9 


12.0 


0.92 


7.2 


43 


9 


534 


0.25 


0.618 


79.2 


1.4 


63.8 


2.9 


5.60 


5.9 


31.7 


5.9 


385 


1.6 


567 


0.4 


0.316 


12.4 


1.5 


8.0 


2.5 


0.72 


8.8 


32.3 


4.3 


172 


1.7 


288 



considered particles) of a collision of a particle with a planet during the lifetime of the particle. 

We define T=Ty,/N as the mean time during which the perihelion distance g of a particle was 
less than the semi-major axis of the planet and Tj as the mean time spent in Jupiter-crossing 
orbits. Below, Psun is the ratio of the number of particles that collided with the Sun to the total 
number of particles. T^*" and T^"^ are the minimum and maximum values of the time until 
the collision of a particle with the Sun, and T^^q and T^q^ are the minimum and maximum 
values of the time when the distance between a particle and the Sun reached 2000 AU. With 
/3>0.01 the values of Pr=l(f'P, T, Tj, T^"'", T^""^', T^^^^, and T^^g are shown in Table 1 for 
A?^=250 (for /?=0.1 we present two runs with 250 different particles), and Psun was obtained for 
all considered particles at a fixed (3. 

The minimum time T™"*" needed to reach the Sun is smaller for smaller particles, but the ratio 
rpmax^j^min jg ^^^^i greater and the ratio Ta^oo/^s'""' is smaller for /3>0.25 than for /3<0.1. For 
/3=0.05, 498 particles collided with the Sun in less than 0.089 Myr, but two particles (with initial 
orbits close to those of the asteroids 361 and 499), which reached 2000 AU, lived for 0.21 Myr 
and 19.06 Myr, respectively. The latter object's perihelion was near Saturn's orbit for a long 
time. At /3=0.05 the first 250 particles did not migrate outside Jupiter's orbit, so Tj=0 in Table 
1. 

For smaller particles (i.e., those with larger (3), Pgun is smaller and the probability of collisions 
of particles with the terrestrial planets is smaller. The probability of a collision of a migrating 
dust particle with the Earth for /3=0.004 is greater by a factor of 1420 than for /3=0.4. These 
probabilities of collisions are in accordance with cratering records in lunar material and on the 
panels of the Long Duration Exposure Facility, which showed that the mass distribution of dust 
particles encountering Earth peaks at d=200 /Ltm [9]. 

ORBITAL EVOLUTION OF ASTEROIDAL PARTICLES 

Several plots of the distribution of migrating asteroidal particles in their orbital elements 
and the distribution of particles with their distance R from the Sun and their height h above 
the initial plane of the Earth's orbit are presented in Figs. 1-7 (Fig. 8 present the results for 
kuipcroidal dust particles). For Fig. 1-2 the number of bins in the semi-major axis a is 1000. 
For the other figures the number of bins in a or i? is 100, and the number of bins in e or /i is 
usually slightly less than 100. The width of a bin in i is 1° (0.5° for Fig. 8). 

In Figs. 1-2, 4-8 to calculate the orbital elements we included the coefficient (1-/3) on the mass 
of the Sun, because of radiation pressure. Fig. 3 is the only exception. Although similar to Fig. 
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4, Fig. 3 (see also figures in [20]) shows the osculating orbital elements that are calculated using 
the normal mass of the Sun (i.e., to transfer from rectangular coordinates to orbital elements 
we used the same formulas as those for massive bodies). For transformations between semi- 
major axes and eccentricities in the above two systems of orbital elements we used the formulas 
obtained by Kortenkamp and Dermott [9]. In total, Figs. 3 and 4 are similar, but they differ in 
some details. For example, for /9=0.4 in Fig. 4 there are some pairs of a and e corresponding to 
a perihelion near Saturn's orbit. There are no such pairs in Fig. 3 for /3=0.4. 

We obtained that the mean time ta (the total time divided by the number N of particles) 
during which an astcroidal dust particle had a semi- major axis a in an interval of fixed width is 
greater for smaller [3 at semi-major axes a<3 AU (exclusive of the gap at a=l AU and /3=0.01). 
In Fig. lb curves plotted at 40 AU are (top-to-bottom) for (5 equal to 0.25, 0.4, 0.05, 0.1, and 
0.01. For /?<0.1 the values of ta are much smaller at a>3.5 AU than at l<a<3 AU, and they 
are usually a maximum at a~2.3 AU. For /3=0.01 the local maxima of ta corresponding to the 
6:7, 5:6, 3:4, and 2:3 resonances with the Earth are greater than the maximum at 2.4 AU. There 
are several other local maxima (Fig. Ic) corresponding to the n:(n-|-l) resonances with Earth 
and Venus (e.g., the 7:8 and 4:5 resonances with Venus). The trapping of dust particles in the 
n:(n-f-l) resonances cause Earth's asteroidal ring [21], [22]. The greater the the smaller the 
local maxima corresponding to these resonances. At /3<0.1 there are gaps with a a little smaller 
than the semi-major axes of Venus and Earth that correspond to the 1:1 resonance for each; the 
greater the /3, the smaller the corresponding values of a. A small gap for Mars is seen only at 
/3=0.01. There are also gaps corresponding to the 3:1, 5:2, and 2:1 resonances with Jupiter. 

For all considered /?, ta decreases considerably with a decrease of a at a<l AU and usually 
decreases with an increase of a at a>5 AU (Fig. 1). Relatively large values of ta at a>40 AU 
for /3=0.05 are due to one particle. For a>5 AU the values of ta are usually a little greater 
at /3=0.25 than those at /3=0.4. The number of particles with a>5 AU at /3=0.25 is smaller 
than at /3=0.4, but they move more slowly to 2000 AU than at /?=0.4. Analyzing Fig. la, we 
can conclude that larger particles make up a greater proportion of the dust population in the 
terrestrial zone than they do in the asteroid belt. 

In Fig. 2 we present the results obtained for the model without planets. In this case, migration 
outside 5 AU is smaller than in the model with planets and, of course, there are no peaks and 
gaps caused by planets (see Fig. 1). For the model without planets, the values of T^***^ were 
about the same as those presented in Table 1, but the values of T^"^ sometimes were smaller. 
With iV=250 the values of Psun for /3=0.25 and /3=0.4 of 0.908 and 0.548, respectively, are 
greater than for the model with planets. For /3<0.1 all of the particles collided with the Sun. 

We now return to the model with planets. At a<4 AU the maximum eccentricities for /3>0.25 
were greater than those for /3<0.1 (Figs. 4-5). At /9=0.01 some particles migrated into the 1:1 
resonance with Jupiter. For a>10 AU perihelia were usually near Jupiter's orbit (for /3=0.05 
and /3=0.25 also near Saturn's orbit). In almost all cases, the inclinations i<50°; at a>10 AU 
the maximum i was smaller for smaller (3 (Fig. 6). 

In Fig. 5 we present the distributions of dust particles in a and e, and in a and i for /?=0.1. 
The plots on the left were obtained for initial positions and eccentricities close to those of the 
asteroids with numbers 1-250, and the plots on the right, to those of the asteroids with numbers 
251-500 (Figs. 1-4, 6-7 were obtained using all particles). Overall, the left and the right plots 
are similar, but for the left plots particles spent more time outside 5 AU. At /3=0.05 none of the 
first 250 particles reached Jupiter's orbit, but two particles with numbers 361 and 499 migrated 
to 2000 AU with perihelia near the orbits of Jupiter and Saturn (Fig. 4). 

Usually there are no particles with h/R>0.7 at R<W AU, with h/R>0.25 at i?>20 AU for 
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Fig. 1. Mean time ta (the total time tasum divided by the number N of particles) during which an 
asteroidal dust particle had a semi-major axis in an interval with a width of (a) 0.005 AU, (b) 0.1 
AU, or (c) 0.001 AU. The values of ta at 1.5 AU are greater for smaller (3. Curves plotted at 40 AU 
are (top-to-bottom) for (5 equal to 0.25, 0.4, 0.05, 0.1, and 0.01. For calculations of orbital elements 
we added the coefficient (1-/3) to the mass of the Sun, which is due to the radiation pressure (the 
same orbital elements are used in Figs. 1-2, 4-8). Initial velocities and coordinates of dust particles 
were the same as those of the first A?^=500 numbered asteroids (Ar=250 for /3=0.01, AT^lOO for 
^=0.004) at JDT 2452500.5. 
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Fig. 2. Same as for Fig. 1 but without gravitational influence of planets and for a smaller number 
of initial dust particles {N=250 at /3>0.05, and iV=100 at /3=0.01). The values of a>6.2 AU were 
reached only at /3=0.4. 



P<0.1, and with h/R>0.5 at R>50 AU for /3>0.25. For /3>0.25 at i?<1000 AU the entire region 
with h/R<0.3 was not empty (Fig. 7). 

The total time spent by 250 particles in inner-Earth (Q = a(l + e) < 0.983 AU), Aten (a < 1 
AU, Q > 0.983 AU), Apollo (a > 1 AU, = a(l - e) < 1.017 AU), and Amor (o > 1 AU, 
1.017 < q < 1.300 AU) orbits was 5.6, 1.4, 4.5, and 7.5 Myr at /3=0.01 and 0.09, 0.08, 0.48, 
and 0.76 Myr at /3=0.4, respectively. The spatial density of a dust cloud and its luminosity (as 
seen from outside) were greater for smaller R. For example, depending on /3 they were by a 
factor of 2.5-8 and 7-25 (4 and 12-13 at /3<0.05) greater at 1 AU than at 3 AU for the spatial 
density and luminosity, respectively. This is in accordance with the observations for the inner 
solar system: inversion of zodiacal light observations by the Helios spaccprobe revealed a particle 
density n{R) oc R~^-^, Pioneer 10 observations between the Earth's orbit and the asteroid belt 
yielded n{R) oc R~^'^, and IRAS observations have yielded n{R) oc R~^'^ [23]. The intensity / 
of zodiacal light falls off with heliocentric distance R as I'^R~'^, with 7=2 to 2.5 and beyond 
about 3 AU zodiacal light was no longer observable above the background light [17]. As in [2], 
we approximately defined the brightness of each particle as R~^. Beyond Jupiter's orbit even 
the number of astcroidal particles at some distance R from the Sun is smaller for greater R, so 
asteroidal dust particles cannot explain the constant spatial density of dust particles at R^3-18 
AU. At such distances, many of the dust particles could have come from the trans-Neptunian 
belt or from passing comets. 

MIGRATION OF KUIPEROIDAL DUST PARTICLES 

We also began a series of runs in which the initial positions and velocities of the particles were 
the same as those of the first trans- Neptunian objects (JDT 2452600.5), and our initial data 
were different from those in previous papers. We store orbital elements with a step of 100 yr. 
At the present time these runs have been finished for (3=0.1 and /3=0.2 with N=50. 

The distributions of dust particles with semi-major axis a and eccenticity e or inclination i, 
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Fig. 3. Distribution of asteroidal dust particles with semi-major axis and eccentricity (designations 
of the number of particles in one bin are the same in Figs. 3-7). For the transfer from rectangular 
coordinates to orbital elements we used the same formulas as those for massive bodies. For Figs. 
1-2, 4-8 for calculations of orbital elements we added the coefficient (1-/3) to the mass of the Sun, 
which is due to the radiation pressure. 
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Fig. 4. Distribution of asteroidal dust particles with semi-major axis and eccentricity (designations 
of the number of particles in one bin are the same in Figs. 3-7). 
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Fig. 5. Distribution of asteroidal dust particles with semi-major axis and eccentricity, and with 
semi-major axis and inclination for /3=0.1. The left figures correspond to the runs with initial 
positions and velocities close to those of the first 250 numbered main-belt asteroids, and the right 
figures correspond to the runs with initial positions and velocities close to those of the asteroids with 
numbers 251-500 (JDT 2452500.5). 
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Fig. 6. Distribution of asteroidal dust particles with semi-major axis and inclination. 
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Fig. 7. Distribution of asteroidal dust particles with distance R from the Sun and height h above 
the initial plane of the Earth's orbit (designations of the number of particles in one bin are the same 
in Figs. 3-7). 
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Table 2. Values of T, Tj, T^^"", T^"^, T^^, T^^g (in Kyr), P^, and Psun obtained for kuiperoidal 
dust particles at /3=0.1 and /3=0.2 (Venus=V, Earth=E, Mars=M) 
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and with distance R from the Sun and height h above the initial plane of the Earth's orbit are 
presented in Fig. 8. The left plots were obtained for /?=0.1, and the right plots were made for 
/3=0.2. For both values of /3, particles with e>0.5 had their perihelia mainly near the semi-major 
axis of the giant planets. The mean eccentricity of kuiperoidal dust particles increased with 
a at a>50 AU, and it exceeded 0.5 at a>60 AU. Bodies that migrated inside Neptune's orbit had 
eOT<0.2. The inclinations were usually less than 35°. Similar figures for intermediate stages of 
the runs for kuiperoidal particles were presented in [24] (for /3=0.05, N=50, and a time interval 
of 2 Myr; and for 13=0.1, iV=100, and t<3.5 Myr). 

At /5=0.4 most particles were outside 50 AU after only 0.02 Myr. If we consider positions of 
particles for t<0.05 Myr, then 84% of them had i?>100 AU. At that time the values of ur/R 
(surface density), ur/R^ (spatial density), and ur/R^ (similar to luminosity) were maximum 
at ~40-55 AU. 

In Fig. 8 the values of h were maximum (30 AU) at i?~57 AU for /3=0.1, and were maxi- 
mum (42 AU) at i?~80 AU for /3=0.2. We usually obtained h/R<0.5, and in the zone of the 
Edgeworth-Kuiper belt the ratio of the number of particles at h=kR dropped usually to 10% 
of the number at h=0 at A;=0.1, but sometimes at /c=0.25. 

Table 2 is similar to Table 1, but was obtained for kuiperoidal dust particles. The fraction 
Psun of kuiperoidal particles collided with the Sun was smaller by a factor of 6 than that of 
asteroidal particles for the same /3, but the difference in collision probabilities with Earth and 
Venus for asteroidal and kuiperoidal particles was smaller than 6. This is due to that the mean 
eccentricities of particles near these planets (especially, at /3=0.2) were smaller for kuiperoidal 
particles than those for asteroidal particles. 

The trans-Neptunian belt is considered to be the main source of Jupiter- family comets, which 
can produce much trans-Neptunian dust just inside Jupiter's orbit. Some of these comets can 
reach typical near-Earth objects' orbits [19]. The total mass of comets inside Jupiter's orbit is 
much smaller than the total mass of asteroids, but a comet produces more dust per unit minor 
body mass than an asteroid. 

MIGRATION OF DUST PARTICLES FROM COMET ENCKE 

We also investigated the migration of cometary dust particles for /3 equal to 0.002, 0.004, 
0.01, 0.05, 0.1, 0.2, and 0.4. For silicate particles such values of /3 correspond to diameters equal 
to about 200, 100, 40, 9, 4, 2, and 1 microns, respectively. The initial positions and velocities 
of the particles were the same as those of Comet 2P Encke. We considered particles starting 
near perihelion (runs denoted as Ato = 0), near aphelion {Atg = 0.5), and when the comet had 
orbited for Pa/4 after perihelion passage, where Pa is the period of the comet (such runs are 
denoted as Atp = 0.25). Variations in time r when perihelion was passed was varied with a 
step 0.1 day for series 'S' and with a step 1 day for series 'L'. For each /? we considered iV=101 
particles for "S" runs and 150 particles for "L" runs. 
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Fig. 8. Distribution of trans-Neptunian (kuiperoidal) dust particles with semi-major axis a and 
eccenticity e or inclination i, and with distance R from the Sun and height h above the initial plane 
of the Earth's orbit. 
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Table 3. Values of T, Ti (in Kyr), and Pr (Venus=V, Earth=E, Mars=M) 
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E 


M 










V 


E 


E 


M 




p 






Pr 


Pr 


T 


Pr 


Ti 




Ato 




Pr 


Pr 


T 


Pr 


Ti 


0.002 





S 


470 


200 


94.6 


12 


551 


0.002 





L 


632 


208 


93.6 


14 


370 


0.002 


0.25 


s 


408 


156 


84.3 


13 


20 


0.002 


0.5 


S 


437 


190 


87.2 


13 


240 


0.004 





s 


370 


148 


62.9 


8.9 


213 


0.004 





L 


303 


139 


66.0 


9.0 


164 


0.004 


0.25 


s 


430 


160 


55.0 


9.3 


109 


0.004 


0.5 


S 


235 


140 


56.3 


8.1 


108 


0.01 





s 


191 


105 


25.1 


5.4 


67 


0.01 





L 


386 


163 


28.5 


6.4 


80 


0.01 


0.25 


s 


238 


86 


24.2 


4.2 


59 


0.01 


0.5 


S 


123 


56 


22.6 


3.8 


49 


0.05 





s 


120 


59 


9.3 


5.7 


1070 


0.05 





L 


142 


67 


7.8 


2.9 


86 


0.05 


0.25 


s 


37 


20 


4.6 


1.6 


5 


0.05 


0.5 


S 


96 


37 


6.4 


2.3 


21 


0.1 





L 


23 


9.1 


6.1 


0.6 


229 


0.1 


0.25 


S 


22 


8.6 


2.8 


0.6 


3 


0.1 


0.5 


s 


13 


6.6 


2.7 


0.47 


3 


0.2 





L 


7.4 


3.5 


3.6 


0.27 


119 


0.2 


0.25 


s 


20 


4.5 


1.9 


0.39 


2 


0.2 


0.5 


S 


12 


3.2 


1.6 


0.22 


2 


0.4 


0.25 


s 


24 


4.3 


1.3 


0.32 


2 


0.4 


0.5 


S 


13 


3.5 


1.4 


0.22 


2 



The results obtained are presented in Table 3. Again, T is the mean time during which the 
perihelion distance q was less than the semi-major axis of a given planet. In contrast to the 
asteroidal dust particles, the values of T did not differ much between Venus, Earth, and Mars 
for the cometary dust particles. Below, Pr = lO^P and T; is the longest lifetime of a particle in 
series of runs. For some runs at /?>0.2, all particles starting close to perihelion got hyperbolic 
orbits just after starting from the comet. These runs are not included in the Table. 

Ah particles with P < 0.01 for Ato = or with (3 < 0.2 for Ato > 0.25 collided with the Sun. 
Relatively large values of Ti for 'S' runs aX /3 = 0.05 and Atg = and for 'L' runs at /3 = 0.1 
and Ato = were due to the particles that reached 2000 AU from the Sun. The values of Pr, 
T, and T/ arc greater for larger particles (i.e., for smaller (3). The values of P,- arc greater for 
Venus than for Earth by a factor of 2 or more. Collision probabilities with Earth were greater 
by a factor of 10-20 than those with Mars and greater for particles starting at perihelion than 
aphelion. For the same values of /3, the probability of cometary dust particles colliding with a 
terrestrial planet was several times smaller than for asteroidal dust particles, mainly due to the 
greater eccentricities and inclinations of the cometary particles. This difference is greater for 
larger particles. 

The ratio of total times spent by cometary particles in inner-Earth, Aten, and Apollo orbits 
was about 1.5 : 1 : 2, but varied from run to run. 

CONCLUSIONS 

We investigated collision probabilities of migrating asteroidal, cometary, and kuiperoidal dust 
particles with the terrestrial planets during the lifetimes of these particles. These probabilities 
were considerably greater for larger asteroidal and cometary particles, which is in accordance 
with the analysis of microcraters. The probability of collisions of cometary particles with the 
Earth is smaller than for asteroidal particles, and this difference is greater for larger particles. 
Almost all asteroidal particles with diameter d>4 /ini collided with the Sun. In almost all cases, 
inclinations i<50°, and at a>10 AU the maximum i was less for larger asteroidal particles. 
The spatial density of asteroidal particles decreases considerably at a>3 AU, so the fraction 
of asteroidal particles among other particles beyond Jupiter's orbit is small. The peaks in the 
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distribution of migrating asteroidal dust particles with semi-major axis corresponding to the 
n:(n+l) resonances with Earth and Venus and the gaps associated with the 1:1 resonances with 
these planets are more pronounced for larger particles. 
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